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Abstract 

We investigate the nonequilibrium tube length fluctuations during the relax- 
ation of an initially stretched, entangled polymer chain. The time-dependent 
variance a 2 of the tube length follows in the early-time regime a simple uni- 
versal power law a 2 = Ay/t originating in the diffusive motion of the polymer 
segments. The amplitude A is calculated analytically both from standard 
reptation theory and from an exactly solvable lattice gas model for reptation 
and its dependence on the initial and equilibrium tube length respectively is 
discussed. The non-universality suggests the measurement of the fluctuations 

(e.g. using flourescence microscopy) as a test for reptation models. 
PACS: 83.20.Fk, 05.70.Ln 
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The most successful framework for the theoretical understanding of systems of entangled 
polymers is the reptation theory developed by de Gennes, Edwards and Doi |TJ,0. One 
imagines the large-scale motion of an entangled flexible polymer as being confined to a tube 
which, loosely speaking, is the sequence of pores occupied by the polymer chain within 
the entanglement "network", made up by the surrounding polymer chains. The motion of 
the polymer transverse to the tube direction is strongly suppressed due to the topological 
constraints imposed by the surrounding polymers. However, mass transport along the tube 
is possible by diffusion of "defects", i.e., folded polymer segments stored within the pores. 
The result of these dynamics is a snakelike random motion where polymer segments may 
retract into the existing tube at its ends (hence effectively shortening it) and subsequently 
diffuse into a new neighboring pore (hence stretching the tube and changing the tube contour 
at the end), while the bulk of the tube does not change its shape. Because of the defects the 
fluctuating equilibrium tube length A* is less than the actual polymer length M of a fully 
stretched chain. 

The basic notions of this picture have been verified experimentally, see @ for the mea- 
surement of the entanglement distance 5 implied in the tube concept, and for the role 
of the surrounding chains as topological constraints. Also extensive molecular dynamics 
simulations confirm the reptation model 0. The theoretical status of some dynamical con- 
sequences, however, remains rather unsatisfactory, a well-known example being the scaling 
behavior of the viscosity 77 oc M K for long flexible chains. Doi has suggested that the dis- 
crepancy between the empirical dependence k rs 3.4 and the reptation result k = 3.0 be due 
to equilibrium tube length fluctuations [0, but this explanation has not remained unchal- 
lenged HH]. To this one may add that also the mass dependence of the effective exponent k 
obtained by Doi appears to be rather severe [HJ. Since within the Doi picture such fluctu- 
ations alter also other observable quantities through a modification of relaxation times (in 
particular, the tube renewal time) we wish to address the question to which extent standard 
reptation theory with underlying Rouse dynamics for individual polymer segments provides 
a proper description of tube length fluctuations. Inspired by the fact that with present 
experimental techniques it is possible to observe and measure the nonequilibrium tube re- 
laxation process for single entangled polymer chains M, we go beyond the older equilibrium 
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analysis of Refs. |||] and calculate here the time-dependent variance of the tube length 
under far-from-equilibrium initial conditions mimicking those prepared experimentally. 

The direct observation of the reptation dynamics was achieved in a series of remarkable 
experiments by Perkins et al. who have studied entanglement effects in a dense solution of 
long A-phage DNA strands using flourescence microscopy || . Attaching a polystyrene bead 
to one end of flourescence-marked DNA molecule (entangled in a background of unmarked 
DNA) and pulling the bead with optical tweezers the DNA was brought into a nonequilibrium 
conformation with an elongated contour directly visible under a microscope. Since this 
contour describes the coarse-grained path of the polymer chain it may be identified with the 
tube. After stopping the bead the contour started to contract along its own path (i.e., the 
tube) as predicted by the reptation model. Since the bead was larger than the estimated 
entanglement distance it could not move and the contraction occured only at the opposite 
free end of the DNA chain. Since the DNA chain used in the experiment was very long (about 
160 times its Kuhn length) it may be regarded as sufficiently flexible to use reptation theory 



as an approximation to study its tube dynamics. In UTTJ we showed that within an universal 
initial time regime (but after a very short transient relaxation process) the measured tube 
length contraction is in agreement with the predictions from both standard reptation theory 
and an exactly solvable lattice gas model for reptation, related to the Verdier-Stockmeyer 
model |12| and the repton model of Rubinstein [O]. 



In order to calculate fluctuations of the tube length from reptation theory we recall 
that within the underlying Rouse theory a polymer chain with bond length b is assumed 
to consist of iV = M/b frictional units, called Rouse segments, which are connected by 
harmonic springs with mean square separation b 2 . By not considering the actual position 
of the polymer in the network but focussing only on the contour length A of the tube the 



reptation dynamics may be cast in a Langevin equation [14 



> 9 , , 3kj^T d 2 , , , > 

C^s(x, t) = b2 dx2 s(x, t) + f(x, t) (1) 

for the position s(x,t) of the Rouse segments in the tube. Here ( is the friction constant 
of the Rouse segments and T is the temperature. Randomness is accounted for by a delta- 
correlated Gaussian random force f(x, t) with zero mean and ( f(x, t)f(x', t') ) = 2(kT5(x — 
x')5(t — t'). By definition the mean tube length is given by A(t) = ( s(N, t) — s(0, t) ). In the 



experimental set-up considered above the position of one end is kept fixed. Hence we have 
to solve ([TJ) with boundary conditions s(0, t) = and the entropic tensile force s'(x, t) = f e 
at x = N for all time t > 0. The entanglement distance S enters only indirectly in so far as 
M S> 5 ^> b is assumed. In the absence of noise ([[]) has the equilibrium solution s*(x) = xf e . 
This yields the mean equilibrium tube length A* = Nf e . 

In order to study the tube length dynamics it is convenient to consider s(x, t) = 
s(x,t) — s*(x) and to introduce the inverse time scale w = 3kBT/((b 2 ). Solving first for 
the homogeneous equation with / = and then incorporating the special solution of the 
inhomogemeous equation we obtain the full solution of the initial value problem 



s(x, t) = ^2 e wp " 1 sin p n xx 
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with the Rouse "momenta" p n = n(2n + 1)/(2N), the sine transform u n (0) = 
1/N J dy sm(p n y)s(y,0) of the initial data and the sine transform g n (r) of the noise 
term f(x,r). For homogeneous initial stretching with initial tube length Ao one has 
s(y, 0) = (Aq — A*)y/N and obtains the mean non-equilibrium tube length 



x „-(2n+l) 2 t/(4T fl ) 



A( f ) = A- + (Ao-A-)^E V +1)2 ■ < 3 > 
with the Rouse time = N 2 / (wtt 2 ). At early times t <Ctr one finds a universal power law 



behavior A(t) = A — 2(A — A*)/n 3 ^ 2 ^t/T R The origin of the scaling exponent 1/2 
is the diffusive motion of the polymer segments. The reduced relaxation function Ai(t) = 
(A - A(£))/(A - A*) depends only on time. 

The tube length fluctuations are determined by the noise term in (0). Specifically for 
the variance cr 2 (t) = ( s 2 (N, t) ) — ( s(N, t) ) 2 one finds 

„, N Nb 2 ( 8 ~ e -( 2 " +1 ) 2t /( 2 ^)\ 

In the prefactor one recognizes the well-known reptation result a* 2 = Nb 2 /3 in equilibrium 
0. We point out that the amplitude depends neither on the initial stretching, i.e. on 
A , nor on the equilibrium tube length A*. The term in the brackets yields the growth of 
fluctuations due to the randomness of the time evolution. At early times 
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which - not surprisingly - exhibits the same universal power law dependence as the tube 
length relaxation. Defining in analogy to the reduced relaxation function the quantity 
A2(t) = cr 2 (t)/a* 2 we note the relation A2(t) = Ai(2£) and, for early times, the simple 
relaxation ratio 

A 2 (t)/A 1 (t) = v / 2 (6) 

which does not contain any parameters of the Rouse-based reptation theory. This feature is 
characteristic for universal amplitude ratios not uncommon in scaling theory, but - as shown 
below - here it is an artefact of the reptation model. 

The result (0) for the tube length relaxation has been recovered within a lattice gas 
approach in which the reptation dynamics are described by the center of mass motion of 
L = M/S individual polymer segments (reptons) of unit length S which we take as the 



experimentally accessible mean entanglement distance in the melt or dense solution ||TT 
The continuous diffusive motion within the tube is modelled by a random hopping of these 
segments between consecutive pores with an exponential waiting time distribution with mean 
Tq- This elementary time scale is the mean first passage time of the motion of one repton to 
a neighboring pore of the tube. To describe the interaction of the polymer with the network 
and with itself we impose the constraint that the consecutively labeled reptons may not pass 
each other within the tube. End reptons may move freely to new pores but are of course 
not allowed to detach from the polymer. 

We assign the label 'particle' to a repton that connects two neighboring pores and the 
label 'hole' to a repton that is fully contained in a pore, i.e., to a "defect" in the language of 
de Gennes (Fig. 1). The diffusion of the defects then corresponds to particle-hole exchange 
with rate l/(2r ), i.e. to particle hopping along the chain with hard-core on-site repulsion. 
The end-point dynamics of the tube correspond to injection and absorption of particles at the 
boundaries. These dynamics are equivalent to the well-known one-dimensional symmetric 



exclusion process with open boundaries [15|. In contrast to Verdier/Stockmeyer and 
Rubinstein resp. [13j we treat one boundary of the particle system (representing the free end 
of the polymer chain) as connected to a reservoir of fixed density p* = A*/M from which 



particles are injected with rate a = p* j (2r ) and absorbed with rate 7 = (1 — p*)/ (2r ). The 
fixed end is described by a reflecting boundary with no particle exchange. We stress that 
the treatment of noise in this approach is rather different from the assumptions inherent in 
the Langevin description ([!]). In particular, because of particle interactions in the lattice 
gas the noise depends of the local hole ( "defect" ) density and hence implicitly on p* which 
determines the equilibrium particle density. 

In this mapping the mean tube length is given by A = 5(N) where N = J2k n k is the 
total number of particles on the lattice and n& = 0, 1 is the particle number on site k. 
Correspondingly, a 2 = 5 2 ((N 2 ) — (N) 2 ) where averages are taken over many realizations 
of the stochastic dynamics, starting from some given nonequilibrium initial distribution 
of tube lengths. An initially fully stretched polymer chain corresponds to an initially fully 
occupied lattice. The Markov generator of this stochastic lattice gas is known as the quantum 
Hamiltonian of the spin-1/2 ferromagnetic Heisenberg chain with boundary field. This is an 



integrable model which can be fully solved by a dynamical matrix product ansatz , similar 
in spirit to the algebraic Bethe ansatz. The defect diffusion as described by the exclusion 
process is the classical analog of the quantum spin waves in the Heisenberg ferromagnet 
and one obtains expressions for expectation values ( (t) . . . rik r (t) ) in terms of sums over 



Bethe wave functions ||17|| . Since we are interested in the early time regime t <^ tr during 
which the bulk conformation remains unchanged we may treat both boundaries separately. 
The tethered end of the chain (caught up in the entanglement) does not contribute to the 
dynamics and the particle number fluctuations induced by the free end can be calculated 
by taking the thermodynamic limit L — > 00. The Bethe sums then turn into integrals over 
pseudomomenta (the analogues of the Rouse modes). 

The calculation of the local density Pk(t) = (rik(t)) and hence of A(t) reduces to the 
solution of a relatively simple random walk problem (j ij and one finds pk{t) = 1 — (1— — 
Z)m=i ex P (~ t/ T o)Ik-m(t/To)]- Here Ik(x) are the modified Bessel functions. By identifying 
w = 1/(2to) one recovers the Rouse result (^) in the regime r <^t <^ tr. 

The calculation of the variance requires the knowledge of the two-point correlation func- 
tion Ck,i(t) = (nk(t)ni(t)) — (rik(t))(ni(t)) from which one obtains a 2 = 5 2 [J2k Pk(t)(l — 
Pk(t)) + 2J2k J2i>k Ck,i(t)}. Computing the correlation function from the Bethe wave function 
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yields after some manipulations the exact expression 



Ck,l(t) — Fk+l-lfe-lit/To) + Fk+l t k-l+l(t/To) 

-F fc _i_i ifc+i (t/r ) - F k -i tlt+ i +1 (t/To) (7) 



with 



4(1 - p*y r*/* 



X 



F mn {x) = -^ — I d6 cos 6 I dve 2vcos2e x 

7T J JO 

e~ 2x I m (2(x-v) cos 6) cos [v sin (29) +n9]. (8) 

Taking the Laplace transform of this expression, summing over the lattice indices and ex- 
pressing the result in terms of elliptic integrals finally yields the asymptotic law (r <C t <C 

° V) = a ^3/2 \ — ( 9 ) 

with the equilibrium variance cr* 2 = <5 2 Lp*(l — p*) of the lattice gas model and g(p*) = 1 + 
( 1 — p* ) (3 — 2 v2) / p* . The details of the rather lengthy calculation will be presented elsewhere 
||18|| . From @ we recover the expected scaling exponent of the growth of fluctuations in 
time. However, considering a different initial conformation |1| reveals that unlike in the 



standard reptation approach the amplitude depends not only on equilibrium quantities, but 
also on the initial tube length. Furthermore, the relaxation ratio of the fully stretched chain 

A 2 (t)/A 1 (t)=g(p*) (10) 

turns out to be a non- universal function depending on the equilibrium defect density 1— p* = 
1 - A*/M (Fig. 2). 

Our main results are the establishment of the universality of the tube length relaxation 
exponent, the proof of nonuniversality of the relaxation ratio g, and the unexpected lack of 
dependence on initial conditions of the amplitude of the variance u 2 (t) in standard reptation 
theory. While the lattice gas is a phenomenological model as is standard reptation theory, 
it is somewhat more microscopic in nature. Unlike the bead-spring model underlying (|]) 
it incorporates the impossibility of stretching the polymer chain and therefore lacks the 
physically problematic parameter b. Hence we believe the master equation describing the 
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lattice gas dynamics to be more trustworthy than its coarse-grained Langevin counterpart. 
This view is confirmed by the amplitude of the tube length fluctuations which vanish in 
the lattice gas description for vanishing equilibrium defect density 1 — p* as they must for 
a chain without defects (whose motion at the boundaries is responsible for the tube length 
fluctuations). In standard reptation theory the amplitude turns out to be unrelated to the 
defect density which points to a conceptual flaw in the treatment of the fluctuations. We 
conclude that the earlier notion that a significant modification of the equilibrium fluctuation 
theory is required extends to non-equilibrium dynamics, in particular, to the treatment of 
noise. Our lattice gas approach suggests that one should incorporate the entropic tensile 
force in the noise amplitude (via multiplicative noise) in order to capture the effect of chain 
connectivity. 

We remark that the experimental data of Ref. || are not suitable for a quantitative 
comparison with our results. On a qualitative level, however, they lend support to the small 
magnitude of the theoretical expression (|9]). In similar experiments specifically designed for 
the purpose of measuring tube length fluctuations, the model dependence of this quantity 
may serve as a test for the reliability of reptation models. 
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Figure captions: 



Fig. 1: Reptation dynamics in terms of the symmetric exclusion process. The polymer 
conformation in the tube (of diameter d ~ 5) is divided into twelve segments of length 8 
occupying six pores corresponding to the particle configuration (from left 

to right). 

Fig. 2: Exact relaxation ratio g(p*) calculated from the symmetric exclusion process. The 
constant at g — y/2 is the standard reptation result. 
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